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I.  INTRODUCTION 

\i 

Algorithms  for  use  in  the  analysis  of  coded  pulse  (band- 
width-  pulse  length  product  greater  than  unity)  processing  systems 
have  been  developed.  These  algorithms  can  be  used  to  delineate 


I desirable  characteristics  of  the  signals  to  be  used  in  the 

processing  system  and  are  of  particular  utility  in  the  digital 
■ simulation  of  system  operation.  In  this  memorandum  the  algorithms 

I are  described  and  examples  of  their  applications  are  given. 
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II.  ALGORITHM  DEVELOPMENT 

A spectral  representation  of  a finite  duration  time  func- 
tion, in  its  usual  mathematical  sense,  consists  of  constructing 
a description  of  the  given  function  in  terms  of  a continuous 
distribution  of  sine  waves  of  infinite  duration.  The  members 
of  this  distribution  are  selected  to  have  amplitudes  and  phases 
such  that  they  add  to  zero  for  all  times  except  during  the 
interval  in  which  they  represent  the  given  function,  in  accordance 
with  Fourier  transform  theory. 

For  purposes  of  crude  analysis,  it  is  often  useful  to 
use  as  a unit  the  finite  duration  sine  wave  instead  of  the  one 
of  infinite  duration.  It  is  only  necessary  to  remember  that  in 
the  final  interpretation  of  the  results  each  line  component  of 

a spectrum  is  in  fact  the  usual  ^ type  spectrum  corresponding 

to  the  sine  wave  of  finite  duration.  In  order  to  formulate  a 
description  of  a signal  processor,  such  a procedure  is  particularly 
convenient  because  it  leads  very  simply  to  results  which  display 
the  spectral  distribution  of  the  signal  as  it  progresses  through 
the  processor. 

On  the  basis  of  the  sampling  theorem  and  by  using  the 
finite  duration  sine  wave  as  a building  block^,  we  make  the 
following  plausible  assertion:  Any  time  function  whose  band- 

width is  p and  whose  duration  is  t can  be  represented  by  some 
combination  of  Pt  finite  duration  sine  waves.  The  choice  of  the 
amplitudes,  relative  phases,  frequency  interval,  and  duration  of 
the  waves  depends  on  the  particular  signal  one  wishes  to  represent, 
as  will  become  clear  below.  It  is  of  Interest  to  consider 
several  examples:  (a)  a sine  wave  of  duration  t is  represented 

by  a single  sine  wave  of  the  same  amplitude,  frequency,  phase, 
and  duration;  (b)  a pulse  of  random  noise  of  duration  t and 


^Alan  A.  Winder  and  Charles  J.  Loda,  j^troduction  to  Acoustical 
Space- time  Information  Processing,  oWl  keport  ACR-63,  pp  46-5o. 


TRACOR,  INC. 


1701  Guodolup#  Sf 


Au9ftn  I,  Tpxos 


bandwidth  ^ is  represented  by  ^T  sine  waves,  all  present  at  one 

time,  at  frequency  intervals  of  duration  t,  an  amplitude 

distribution  corresponding  to  the  power  spectrum  of  the  pulse, 

and  random  phases;  (c)  an  FM  slide  of  duration  t and  bandwidth 

P is  represented  by  Pt  sine  waves  of  duration  with  one  such 

wave  present  at  any  one  time,  the  amplitude  of  each  wave  the 

same  as  that  of  the  pulse,  and  occurring  in  sequence  at  frequency 

intervals  — in  the  order  prescribed  by  the  direction  of  the 

slide;  and  (d)  a type  of  pseudo-random  noise  pulse  of  constant 

amplitude,  bandwidth  P and  duration  t is  represented  by  Pt  sine 

waves  of  duration  i,  spaced  at  frequency  intervals  of  multiples 
1 ^ 

of  — , with  the  different  frequencies  occurring  in  a random  order 
such  that  only  one  wave  is  present  at  any  one  time.  The  ampli- 
tude of  each  Incremental  wave  is  set  equal  to  the  pulse  amplitude. 

Using  the  above  concepts,  a section  of  random  noise  of 
bandwidth  P and  duration  t is  denoted  by  the  general  form 

M i 

n(t)  = E cos  [2Tr(fQ  + -)  t + 5,^^] 

i=-M 


where  f^  is  the  center  frequency  of  the  band,  M = ^ 

is  an  arbitrary  phase  angle, and  it  is  understood  that  the  sinu- 
soidal waves  all  begin  at  some  time  t = a and  end  at  t = a + t. 

This  form  is  most  useful  for  the  representation  of  CW  and  random 
noise  signals.  Note  that  in  this  form,  the  instantaneous 

amplitude  of  n(t)  may  fluctuate  widely  between  the  limits  of 
M 

t E |n. I , so  that  this  form  is  of  limited  usefulness  in  construct- 
i=-M  ^ 

Ing  a relatively  constant  amplitude  pulse. 

The  second  form  to  be  used  is  better  suited  to  represent 
FM  slides  and  a certain  class  of  pseudo-random  noise.  Here  a 
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section  of  signal  or  noise  of  bandwidth  p and  duration  t is 
represented  by  the  general  form 

a(t)  = I p(i,t)  a.  cos  [2Tt(f,  + t + 

i=l  ^ i T 1 

where  is  the  lowest  frequency  in  the  pass  band,  p(i,t)  is  a 
"pulse"  or  "gating"  function  defined  by 


p(i,t)  = 0 if  pt  < i-1  or  Pt  > i 
= 1 if  i-1  s Pt  < i 

where  m(i)  denotes  m,  an  integer  in  the  set  of  numbers  from  1 
to  Pt  ordered  on  the  index  i,  and  the  phases,  are  chosen 
so  that  the  function  is  continuous.  If  the  signal  is  of  constant 
amplitude,  then  the  phase  angles  are  defined  by  the  recurrence 
relationship 


i+1  “ [“(i)  - m(i+l)]  + 


In  the  case  of  an  FM  slide,  the  numbers  m(i)  are  chosen  in  order 
with  respect  to  the  index  i,  and  all  the  amplitudes  a^^  are 
equal  to  the  amplitude  of  the  pulse.  For  the  case  of  pseudo- 
random noise,  the  numbers  m(i)  are  chosen  randomly  with  respect 
to  the  index  i and  the  amplitudes  a^  may  be  chosen  arbitrarily. 

Note  that  the  maximum  instantaneous  amplitude  of  the 
pulse  is  given  by  the  maximum  value  of  a^^.  Here  the  peak  ampli- 
tude can  be  easily  controlled,  and  therefore  this  form  is  useful 
in  constructing  a pseudo-random  noise  pulse  of  constant  amplitude 
by  choosing  the  a^'s  to  be  equal. 


' rsfi 


4 


'>^J'  -4iyB|il!!^^.- 


r^i  ->iCByiin<ig^ 
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The  differences  between  the  two  algorithms  are  very 
significant  if  they  are  to  be  used  for  the  design  of  sonar 
pulses.  By  using  the  first  representation  in  which  Pt  sine 
waves  of  duration  t are  added  simultaneously,  the  power  spectral 
density  of  the  generated  signal  can  be  closely  controlled,  while 
the  amplitude  of  the  signal  envelope  may  fluctuate  between  wide 
limits.  By  using  the  second  representation  in  which  Pt  sine 
waves  of  duration  are  combined  sequentially,  the  amplitude  of 
the  signal  envelope  may  be  closely  controlled,  while  the  power 
spectral  density  may  contain  wide  fluctuations. 


The  difference  in  ease  of  control  of  the  envelope  behavior 
as  a function  of  time  can  be  noted  by  examining  the  two  formulariza- 
tions.  In  the  case  of  the  simultaneous  sine  waves,  the  amplitude 
at  any  time  is  determined  by  the  summation  of  Pt  terms,  each 
term  having  an  amplitude  coefficient  modified  by  a phase  angle. 

In  the  case  of  the  sequential  representation,  the  maximum 
amplitude  during  any  period  of  time  is  determined  by  only  one 
coefficient. 


The  difference  in  frequency  behavior  of  the  power  spectral 
density  can  be  noted  by  considering  the  voltage  spectra  of  the 
incremental  sine  waves.  In  the  case  of  the  simultaneous  sine 
waves  of  duration  t , the  spectra  are  quite  narrow  so  that  when 
the  spectrum  of  the  composite  signal  is  obtained  by  adding  the 
individual  spectra  of  each  incremental  wave  the  phases  of 
adjacent  spectra  do  not  interfere.  In  the  case  of  the  sequential 
sine  waves  of  duration  the  spectra  of  the  incremental  sine 
waves  are  quite  wide  and  overlap,  so  that  the  phases  must  be 
considered  in  obtaining  the  composite  spectrum. 

The  remainder  of  this  memorandum  will  be  devoted  to  the 
application  of  these  algorithms  in  working  example  problems. 
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III.  HETERODYNE  CORRELATION  WITH  RANDOM  NOISE 

The  first  representation  outlined  above  can  be  used  to 
indicate  the  performance  of  a heterodyne  correlator.  Let  the 
transmitted  pulse  of  pseudo-random  noise  be  denoted  by  s(t), 
where 

M . 

s(t)  = z cos  [2TT(fo  +“)t  + t-]. 

i=-M 

The  replica,  r(t),  of  the  transmitted  pulse  for  use  in  the  signal 

channel  of  the  correlator  is  then,  after  inverting  the  spectrum  i 

of  s(t),  about  f , 
o 

M . I 

r(t)  = z cos  [2TT(f  - i)  t - iii.].  | 

1=M  u T 1 ^ 

I 

! 

It  is  convenient  to  first  examine  the  action  of  the  i 

correlator  with  an  input  of  random  noise.  Denote  the  received 
noise  by 

M 

n(t)  = E n.  cos[2iT(f„  + -)t+o.  ]. 

i_M  ^ O T 

where  the  received  noise  is  not  correlated  with  the  transmitted 
signal.  Forming  the  product  r(t)  n(t)  = G(t)  yields 

M . M 

' J-M  ""i  -*iJ  "k  <=°n2^(v  + »ki- 

or, 

G(t) 


6 


M 

V. 


M 

E 


i=-M  k=-M 


cos[2TT(2f^  + is: 


Ot  + + 


J 
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M M r.n,  « , . 

+ E E cos[2TT(^)t  + ^ + n ■] 

i=-M  k=-M  T K i 


The  first  double  sum  term  of  G(t)  contains  all  the  sum  frequencies 
of  the  components  of  r(t)  and  n( t) , while  the  second  double  sum 
term  contains  all  the  difference  frequencies.  It  is  convenient 
to  consider  only  the  sum  frequencies  and  to  collect  terms 
corresponding  to  the  same  frequency  by  altering  the  summation 
indices.  Collecting  L'ermI  we  have 

G (t)  =5;  s . cos[2n(2f  +^)t  + 

® j=-2M  i=-j-M  ° ^ ^ 


+ E 
i=-M 


COS[2TT(2f^)  t-hp^  - til 


2M  M-j  • 

cos[2„(2f^  +-l)t  + 

J = 1 X=-M 


By  expanding  the  cosine  in  each  of  the  sums  in  the  above  expression 
it  can  be  rewritten  as 


rr  M r.n.  . _ 

.(t)  = E I E — cos(cp.^.  - i(i.)| 


M r.n 


cos[2n(2f^  +1)  t + q.]  +(  r 2 cos(cp.-^jl 
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P.1 

r ^ 

[ ^ sin(cp^-i|i^)J  J cos[2Tr(2f^j)t  + 


2M  r r M- j r.n.  , -, 


rM- j 4 1 T ' 


cos[2Tr(2f^  + 4)  + 6j3 

where  the  are  constant  phase  angles  for  the  particular  sample 
of  length  T under  consideration.  The  above  function,  Gg(t) 
represents  the  input  to  the  heterodyne  correlator  band-pass 
filter  located  at  2f^(or  the  bank  of  doppler  filters  in  the  case 
of  the  SQS-26(XN-2)) . 

The  following  properties  can  be  deduced  from  the  above 
equation,  where  the  list  also  Includes  those  which  were  obvious 
from  the  beginning  and  have  been  verified  by  more  sophisticated 
analysis . 

1.  The  component  of  Gg(t)  located  at  2f^  consists  of  the 
summation  of  (2M+1)  or  Pt  terms. 

2.  The  summations  consist  of  successively  fewer  terms  as 
the  frequency  gets  further  away  from  2f„.  Thus  for 
frequency  (2f^t  — ) it  consists  of  2M  terms,  for  fre- 

p 

quency  (2f^t  — ) of  (2M-1)  terms,  etc.,  until  for  fre- 

PM 

quency  (2f^t  ~)  it  consists  of  only  one  term. 

3.  All  useful  information  at  this  point  is  contained  in 
the  amplitude  of  each  frequency  component. 
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4.  The  amplitude  of  each  of  the  components  in  the  spectrum 
i.e.,  the  value  of  the  sum  on  i for  given  j,  depends 
on  the  values  of  the  epoch  angles,  completely  analogous 
to  the  situation  encountered  in  beam  forming.  If 
2M  is  reasonably  large  and  the  phase  differences 
^^i+j  " '*'i^  random  as  a function  of  i,  the  value 

of  each  sum  on  the  index  i is,  on  the  average,  for  all 
possible  combinations  of  epoch  angles,  proportional  to 
the  square  root  of  the  number  of  components,  i.e.,  the 
terms  sum  incoherently. 

Now  consider  the  input  of  the  correlator  to  contain  both  a 
zero  doppler  signal  and  the  uncorrelated  noise.  The  received 
waveform  is  then 


n(t)+s(t) 


M 

E 

i=-M 


cos[2n(fQ  + -^)t  + cpj^]+s^  cos[2n(fQ 


where  the  subscript  "s"  on  the  phase  of  the  signal  portion  is  used 
to  note  the  effects  of  possible  changes  between  the  transmitted 
and  received  signal.  Forming  the  product  function  G(t)  = r(t) 
[n(t)  + s(t)j  gives 

M i ^ k 

G(t)  = ^ r.  cos[2Tr(f  - -)t  - i n,  cos[2n(  f _-t-)  t + cpu] 

^ ix  W k 

+ I r.  cos[2TT(f  - -)t:  - >|i.]  z s,  cosr2Tr(f  +-)+  ] 

i__M  ^ o T 1 k=-M  o T KS 

or 
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M M rr.n, 

G(t)  = y T \ — 5-  cos[2Ti(f  + t + cpk-'I'i] 

i=-M  k=-M  I o T K 1 

r s • 

+ ^ cos[2TT(2fo  + ^)t  + - i,.]  j 

M M rr.n,  . . 

M M cos[2n(i^)t+cpi^  + *.] 

i=-M  k=-M  ^ 

IT  S • ^ 

■*■  cos[2n(-^^) t+\|ij^g  - 


where  the  first  double  sum  contains  all  the  sum  frequencies  and 
the  second  double  sum  contains  all  the  difference  frequencies. 
It  is  again  convenient  to  consider  only  the  sum  frequencies  and 
to  collect  terms  corresponding  to  the  same  frequency. 

Collecting  terms,  we  have  for  the  sum  frequency 


-1  M rr.n.  . 

(t)  = y y \ ^ cos[2TT(2f  + J-)t  + . - !);.] 

j=-2M  i=-j-M  o T i+j  1 


+ cos[2Tr(2f^  + ^)t  + - ^i]} 


M r.s. 

cos[2n(2f^)t  + cos[2Tr(2f^)t  + 


r . s . 


2M  M-j  rr.n.  . . 

jh  iLm  1“^  cos[2„(2£„  + i)t  + 


^ • S • • 

+ cos[2TT(2f^  + 'i)t  + - till  . 
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The  component  located  at  2f^  is  of  particular  interest 
because  it  is  the  frequency  which  excites  the  heterodyne  correlator 
band-pass  filter.  Denoting  this  term  by  G^Ct),  we  have 

M rr.n.  r^s. 

Gf(t)  = ^ I 1 cos[4TTfpt  + ^ cos[4TTfQt+  ^is"’''i^r 

i=-M  1 J 


Rearranging  we  have 


2'l| 


|-  r M r .n . 1 ^ r ^ ^i"l  1 

l[i=?M  -r-  sin(«>i-'fi)J  / 


cos[4nfQt+o^] 
M r.s. 


rr  ri  LjO. 

I cos(ili.  - l||  . ) 

Ui=-M  IS  xj 


r ^ 11^ 

E — ^ ®^’^('^is  “ "fi^  I cos[4Trfot  + a .] 

- i=-M  d 


where  0^  and  are  constant  phase  angles  for  the  particular 
sample  of  noise  of  length  t under  consideration. 

Examining  the  above  equation  yields  an  indication  of  the 
processing  gain  of  the  heterodjme  correlator.  The  output  component 
at  2f^  due  to  the  noise  sums  incoherently  because  of  the  random 
nature  of  the  phase  difference  and  thus  the  output  due  to 

noise  is  approximately  proportional  to  the  square  root  of  the  number 
of  terms  (or>/^  since  there  are  Pt  terms).  The  output  component 
due  the  signal  sums  coherently  because  the  phase  difference  ( iji . -f  .) 
is  approximately  zero.  Hence  the  signal  output  is  proportional  to 
the  number  of  terms  (Pt).  Thus  the  amplitude  of  the  component  at 
frequency  2f^  due  to  a signal  is  greater,  by  a factor  proportional 
to  the  square  root  of  the  pulse  duration-bandwidth  product,  than 
when  noise  alone  is  present  at  the  same  power  level.  The  power 
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density  spectra  associated  with  the  above  time  functions  are 
shown  in  Figure  1.  Note  that  because  of  the  finite  pulse  length, 
each  of  the  lines  actually  represents  a sin  x/x  type  spectrum 
so  that  the  composite  spectrum  over  the  band  is  relatively  smooth 


12 


The  algorithm  may  also  be  used  to  determine  the  degradation 
in  performance  introduced  by  using  phase-shift  beamforming  instead 
of  the  correct  time-delay  beamforming.  The  example  transducer 
geometry  is  shown  in  Figure  2.  The  transducer  consists  of  24 
vertical  elements  spaced  at  5®  intervals  along  an  arc  of  an  8 ft. 
radius  circle. 


For  the  case  of  an  incident  wave  arriving  in  the  direction 
as  shown  in  Figure  2,  one  method  of  time  delay  beamforming  would 
be  to  delay  the  output  from  each  element  to  obtain  equivalent 
time  arrival  at  the  line  joining  elements  12  and  12'  (normal  to 
the  direction  of  the  incident  wave).  The  appropriate  time  delay 
for  the  e'th  element  is  given  by 

*"e  “ 566(5  " •5)5°  - cos  57.5°]. 

Summing  the  signals  over  all  of  the  elements,  the  composite  beam 
signal  is  given  by 
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M 12 


s(t)  = 2 r S^e  cos[2Ti(f  +4) 

i=-M  e-1  ° ^ e 1 


i=“M  e=i 


In  phase  shift  beamforming,  the  quantity  (t^f^  + — ^)  is  replaced 
by  its  narrow  band  approximation,  If  it  is  assumed  that 

the  phase  shift  networks  are  designed  such  that  signals  of  frequency 
f^  from  the  different  elements  are  added  in  phase,  then  the  signals 
of  adjacent  frequencies  will  be  in  phase  error.  The  correctly 
formed  signal  from  the  array  will  be  given  by 

M 12  . 

s(t)  = 2 ^ cos[2Tr(fQ+  -^)t  + 

while  the  phase-shift  formed  signal  will  be  given  by 


^ i=?M  e-1  + -)t  + + 2nt^  i]. 

For  simplification,  assume  that  all  elements  have  the  same 
response  in  the  direction  of  the  incident  beam  and  that  the  fre- 
quency spectrum  of  the  signal  is  flat  over  the  bandwidth.  There- 
fore the  s^g  are  all  equal  and  can  be  set  equal  to  unity.  Further 
note  that  the  phase  error  is  doubled  for  two  way  transmission  through 
the  system.  Then  we  have  for  the  correctly  formed  and  the  phase- 
shift  formed  signal  respectively 

s(t)  =24  z cos[2Ti(f  + -i)t  + 

i=-M  o T 1 

and 
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M 12  , i 

®2,('=)  - ® c°s[2n(f„  + -)t  + . 

Rearranging  the  phase-shift  formed  signal  yields 

M r r 12  . r 12  . 

® iLm  I [ eSl  t’J  / 

cos[2n(f^  + l)t  + + yj 

where  j^2 

E sin(4TTtg  -^) 

tan  6.  = 

^12  . 

E cos(4TTt  -) 
e=l  ® ^ 


In  order  to  actually  determine  the  effect  of  the  phase- 
formed  beam  on  system  operation,  consider  the  processing  system 
to  be  a heterodyne  correlator.  Multiplying  the  phase- shift  formed 
signal  with  the  spectrum  inverted  replica  of  the  exact  signal 
yields 


E COS[2TT(f  - ~)t 
i=-M  ° ^ 


12 

E 

e=l 


sin(4ntg  7)]  } 


cos[2n(fQ  + 7)1  + + 8j^] 


Expanding  Gg^Ct),  and  considering  only  the  sum  frequency  located 
3b  2f^,  yields  the  signal  as  would  be  seen  by  the  correlator  band- 
pass filter  at  perfect  registration  between  the  received  signal 
and  the  replica.  Then  this  signal  denoted  by  G£(t),  be  comes 
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Gf(t)  «=  cos(4ntg  i)J  +1"  sin(4ntg  ^)j  | cos(4nfQt  + 0^] 

Because  reduce  the  above  equation  to 


M X2 

Gf(t)  «j^  1 + 2 cos(4ntg  -^)j  + [ sin(4ntg  -^)j  | cos 


cos  4jTf  t. 
o 


Evaluating  the  above  expression  for  the  case  of  a 100  cps 
bandwidth  and  a 0.5  second  pulse  and  comparing  with  the  similar 
results  obtained  using  the  correctly  formed  beam  shows  a performance 
degradation  of  0.2  db  due  to  the  use  of  the  phase-shift  beam  forming. 


t 
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t 
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V.  LINEAR  FM  SLIDE  SIGNAL 

For  Che  case  of  an  FM  slide  of  duration  t and  bandwidth 
P,  the  sequential  frequency  representation  is  preferred.  The 
slide  is  represented  by  Pt  sine  waves  of  duration  ^ with  only 
one  such  wave  being  present  at  one  time,  and  occurring  in 
sequence  at  frequency  intervals  of  in  the  order  prescribed  by 
the  direction  of  the  slide.  The  amplitude  of  each  incremental 
wave  is  the  same  as  the  amplitude  of  the  pulse , and  the  phase  of 
each  incremental  wave  is  chosen  so  that  the  function  is  continuous. 
Thus  an  upward  sliding  transmitted  signal  can  be  represented  by 

P T ' 

S(t)  = a z p(i,t)  cos[2Tr(f^  + ^^)t  + i^.]. 

i=l  1 T 1 j 

i 

While  the  above  representation  is  general,  it  is  convenient  to  ’ 

simplify  it  for  the  specific  case  of  an  FM  slide.  The  upward  j 

sliding  transmitted  signal  can  also  be  represented  by  i 

S(t)  = a cos[2n(fi  + + 'Kk(t)] 

I 

where  the  index  k(t)  is  defined  by  | 

j 

k(t)  s largest  integer  less  than  or  equal  to  Pt,  the  phase 
angle  defined  by  | 

1 

♦k(t)  = i 

J 

j 

and  f^  is  the  initial  frequency  of  the  slide.  i 

The  replica  of  the  transmitted  pulse  for  use  in  the  reference  j 
channel  of  the  heterodyne  correlator  is  obtained  by  inverting  the  I 
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spectrum  of  the  transmitted  pulse  about  the  center  frequency  of 
the  slide.  The  replica  then  becomes 

r(t)  = a cos[2TT(fj^  + P - - *k(t)^* 

Generally  only  a portion  of  the  returned  signal  will  be  In  the 
signal  channel  of  the  correlator.  Thus  It  Is  Instructive  to 
examine  the  behavior  of  the  correlator  output  as  the  received 
signal  slides  Into  the  correlator.  As  the  signal  slides  Into  the 
correlator,  It  Is  present  only  for  duration  a where  o < t.  The 
received  signal  (without  noise)  Is  represented  by 

s(t)  = b(t)  cos[2n(fj^  + ^U^T))(t+g-T)  + 

where  K(t+(j-T)  s largest  Integer  less  than  or  equal  to  (t4tj-T)P, 

b,  r ~ o ^ bs  T 

b(t)  = 

0>  t < T - CT  , 

and  time  Is  measured  from  the  beginning  of  the  replica  pulse  used 
In  the  correlator. 

The  product  function  G(t)  = r(t)  s(t)  yields 
G(t)  = a b(t)  cos|2n[fj^  + ( t4r7  - t ) + 

cos{a<fi  + e - . 
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Expanding  the  product  of  the  cosines  yields 
G(t)  = cos|2TTr2fj^  + P + 

+ 2TT[fj  + - *k(t)} 

+ cos|2tt[-P 

+ 27T[fj^  + + tk(t-hj-T)  + *k(t)}  • 


It  Is  asserted  that  the  first  term  of  the  above  equation 
approximates  a sinusoid  of  amplitude  and  frequency  2fj^  + 

The  validity  of  this  assertion  can  readily  be  demonstrated  at  the 
points  where  t - a = ^ and  b = with  m and  n being  integers. 

Here  k(t)  = Pt  and  k(t+tj-T)  = p(t+g-T).  When  substituted  into  the 
above  equation  the  first  term  (or  the  sum  frequency)  becomes 


GgCt)  = cos  |2TT[2fj^  + p t 

+ 2TT[f^  + P](o-t)  + ^ (Pt+D-  ^ (t+a-T)[P(t4n-T)+l]j 


which  reduces  to 

Gg(t)  = cos[2TT(2fj^+  ^)t  + 2n(fj^  + ^2^f^)(a-T)] . 


The  following  conclusions  can  be  drawn  from  the  preceding  represen- 
tation of  the  FM  slide. 
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If  there  is  no  doppler  shift  and  the  received  signal 
is  in  time  registry  with  the  stored  replica,  the  correla 
tor  output  will  excite  the  zero  doppler  filter  located 
at  a frequency  of  2fj^  + 0. 

If  there  is  no  doppler  shift  and  the  received  signal 
is  not  in  time  registry  with  the  stored  replica,  the 
correlator  output  will  excite  the  doppler  filter  located 
at  a frequency  of  2fj^  + 

If  there  is  a doppler  shift,  the  correlator  output  will 
excite  the  zero  doppler  filter  for  some  value  of  the 
time  registry  a;  hence  only  one  doppler  filter  is 
required  to  determine  the  presence  of  a received  signal. 
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VI.  UNIFORM  AMPLITUDE  PSEUDO-RANDOM  NOISE  I 

2 ^ 

Studies  by  Stewart  and  Westerfield  indicate  that  a noise- 
like signal  may  be  superior  to  either  FM  slides  or  CW  sigials  in  ^ 

several  applications.  The  use  of  an  actual  random  noise  signal  is 
not  particularly  desirable  because  most  sonar  systems  are  peak- 
power  limited  and  random  noise  has  large  amplitude  variations. 

While  it  is  relatively  easy  to  obtain  noise  of  the  desired  bandwidth, 
it  has  been  difficult  to  generate  a noise-like  signal  of  constant 
amplitude  which  also  has  a uniform  energy  density  spectrum  over 
each  relatively  short  duration  signal.  The  sequential  frequency 
representation  (of  which  the  preceding  representation  cf  the  FM 
slide  is  a special  case)  is  well  suited  for  representing  a certain 
class  of  pseudo-random  noise  which  has  desirable  characteristics. 

a 

Consider  the  signal  to  be  composed  of  incremental  sine 

waves,  each  of  duration  i,  occurring  at  frequencies  spaced  at 
1 ^ 

intervals  — over  the  bandwidth,  and  with  only  one  such  incremental 

T 

wave  present  at  any  one  time.  Note  that  the  order  in  which  the 
individual  frequencies  occur  has  not  been  specified.  Thus  a 
transmitted  signal  of  constant  amplitude  can  be  represented  by 


( • \ 

S(t)  = a T.  p(i,t)  cosr2TT(f,  + + \|j . ] 

i=l  i T 1 


where  p(i,t)  is  a "pulse"  or  "gating"  function  defined  by 

0 if  pt  < i - 1 
p(i,t)  = 1 if  i-1  s pt  < i 

0 if  i < pt 


■J.L.  Stewart  and  E.C.  Westerfield,  "A  Theory  of  Active  Sonar 
Detection,"  Proceedings  of  the  I.R.E. , pp  872-881;  May,  1959 • 
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m(i)  denotes  m,  an  integer  in  the  set  of  numbers  from  1 to  Pt 
ordered  on  the  index  i,  and  is  defined  by  the  recurrence  re- 
lationship 

'‘'i+l  " [®(i)  - ni(i+l)]+  ,1,^ 

so  that  the  function  s(t)  is  continuous. 

The  replica  of  s(t)  is  formed  by  inverting  frequency  and 
phase  about  f^^,  so 

r(t)  = S p(i,t)  cos[2TT(fi  - t - iti.l. 

i=l  1 T 1 

Denote  the  received  signal  as  it  slides  into  the  correlator  by 
S(t)  = E p(i,t)  COS[2TT(fi  + t-iT7-T)  + A.  .] 

i=j 

where  the  received  signal  overlaps  the  replica  for  time  a and 
j = P(t  - a). 

Forming  the  product  G(t)  ~ s(t)r(t)  yields 

G(t)  = E p(i,t)  cos[2Ti(f,  -4r^)t  - A-] 
i=j 

COS[2TT(fj^  + -^^~j^)(t+a-T)  + 


Expanding  the  product  of  the  cosines,  and  considering  only  the  sum 
frequency  terms  yields  the  input,  G_(t),  to  the  doppler  filter 

O 

bank.  Hence 
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Gs(t)  = 


Z p(i,t)  cos[2Tr(2f  + - ij, 

i=j  O T 1 J 1 


+ 2TT(f^  + -^~^)(a-T)]. 


Considering  the  range-range  rate  ambiguity  and  examining 
the  above  equation  indicates  the  most  desirable  type  of  signal 
which  can  be  generated  using  this  representation.  Consider  the 
following  types  of  signal  generated  by  the  above  procedure. 

(1)  Assume  that  all  of  the  m's  are  equal.  This  reduces  to 
the  case  of  a CW  signal  which  has  good  range  resolution  but  poor 
frequency  resolution  if  t is  short,  and  poor  range  resolution  but 
good  frequency  resolution  if  t is  long. 

(2)  Assume  that  the  m's  are  chosen  in  a linear  sequence. 
This  reduces  to  the  case  of  an  FM  slide  which  possesses  a range- 
range  rate  ambiguity  in  which  range  can  be  determined  accurately 
only  if  the  range  rate  is  known  and  range  rate  can  be  determined 
accurately  only  if  the  range  is  known. 

(3)  Assume  that  the  m's  are  chosen  randomly.  This  causes 

the  frequency  shift  to  be  a random  function  spread 

over  a range  of  2P  unless  j = 0 when  there  will  be  no  frequency 
shift.  Hence  the  signal  is  capable  of  good  resolution  in  both 
range  and  range  rate.  As  previously  noted,  the  signal  is  very 
desirable  because  it  retains  the  constant  amplitude  property. 

The  algorithm  described  in  this  section  has  been  used  to 
generate  several  constant  amplitude  PRN  signals.  The  characteris- 
tics of  these  signals  were  investigated  by  measuring  their  power 
density  spectra  and  by  measuring  various  correlation  functions. 
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The  measured  power  density  spectra  of  three  PRN  signals  are 
shown  in  Figures  3-5.  In  each  case,  the  signal  had  constant 
amplitude  over  its  duration  of  1/2  second  and  was  intended  to  have 
a flat  frequency  spectrum  over  the  range  of  100  cps  to  200  cps. 

While  the  signals  did  not  have  the  desired  rectangular  power 
spectra,  they  did  have  the  essential  characteristics  for  use  as 
sonar  signals.  Very  little  signal  power  was  outside  of  the 
desired  bandwidth  and  the  power  was  distributed  over  the  entire 
bandwidth,  even  though  the  distribution  was  not  as  uniform  as  could 
be  desired.  For  comparison,  the  power  density  spectra  of  two 
linear  FM  slides  are  shown.  The  spectrum  shown  in  Figure  6 was 
obtained  for  a slide  which  was  generated  by  using  the  algorithm 
of  this  section,  while  the  slide  used  to  obtain  Figure  7 was  generated 
by  using  the  true  mathematical  representation  of  a linear  slide. 

The  correlation  functions  of  these  pseudo-random  noise 
signals  are  shown  in  Figures  8-10.  The  set  of  correlation 
functions  shown  in  each  figure  actually  represent  sections  sliced 
from  an  ambiguity  diagram  with  all  of  the  slices  being  parallel 
to  the  time  axis . The  nine  different  bands  of  each  figure  represent 
the  polarity  cross-correlation  functions  of  the  original  signal 
against  a replica  of  the  signal  which  has  been  modified  by  shifting 
i t by  a constant  frequency  increment.  The  nine  frequency  increments 
used  are  8,  6,  4,  2,  0,  -2,  -4,  -6,  and  -8  cps.  The  center  band 
represents  zero  frequency  shift  and  thus  is  the  auto-correlation 
function . 
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FREQUENCY 

MEASURED  POWER  DENSITY  SPECTRUM  OF  TRUE  FM  SLIDE 
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